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ABSTRACT 

We present new multiwavelength submillimeter continuum measurements of the cir- 
cumstellar dust around 48 young stars in the p Ophiuchus dark clouds. Supplemented 
with previous 1.3 mm observations of an additional 99 objects from the literature, the 
statistical distributions of disk masses and submillimeter colors are calculated and com- 
pared to those in the Taurus-Auriga region. These basic submillimeter properties of 
young stellar objects in both environments are shown to be essentially identical. As 
with their Taurus counterparts, the p Oph circumstellar dust properties are shown to 
evolve along an empirical evolution sequence based on the infrared spectral energy distri- 
bution. The combined p Oph and Taurus Class II samples (173 sources) are used to set 
benchmark values for basic outer disk characteristics: Md ~ 0.005 M©, Md/M^^ ~ 1%, 
and a ~ 2 (where Fi, oc between 350 pm and 1.3 mm). The precision of these num- 
bers are addressed in the context of substantial solid particle growth in the earliest 
stages of the planet formation process. There is some circumstantial evidence that disk 
masses inferred from submillimeter emission may be under-estimated by up to an order 
of magnitude. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — 
solar system: formation — stars: pre-main— sequence 



Introduction 



Circumstellar accretion disks are an integral part of the star and planet formati on process. 
They are critical for conserving angular moment um during molecular cloud co re collapse (jTerebev et al 
1984 ). channeling mass onto the centr al star (ILvnden-Bell &: Pringld 11974 ). and assembling raw 
material into a planetary system (e.g.. iLissauer &: StevensonI 120071 ) . The structural, material, and 
chemical compositions of disks provide crucial insights on the initial conditions of the planet for- 
mation process. While molecular gas is the dominant mass species in these disks, observations of 
this phase are notoriously difficult due to the radiative inefficiency of H2 , the depletion of molecules 
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from the gas phase onto dust grain surfaces, and the relative insensitivity of many current tele- 
scopes. Consequently, much of our knowledge of disk properties is gathered from the comparatively 
bright thermal co ntinuum emission of dust particles heated by irradiation from the central star (see 
Natta et alJboOTl l. 



At any given frequency, this continuum emissi on is generated from d ust grains at a relatively 
small range of contributing stellocentric radii (e.g.. iBeckwith et al.lll99Cll ). Therefore, the detailed 
morphology of its spectral energy distribution (SED) is dictated by the physical structure and ma- 
terial content of the disk (e.g., temperatures, densities, geometry, and particle properties). Much of 
the infrared SED is generated in the optically thick inner few AU of the disk, corresponding roughly 
to the region in the Solar System occupied by the terrestrial planets. Because optic al depths are 



so high at these radii, the infrared SED can be used t o infer disk temperatures (e.g., lAdams et al 



1987 



Chiang &: Goldreichlll997l : iD'Alessio et al.lll998l ). The far-infrared and millimeter /radio por 



tion of the SED (hereafter referred to as "submillimeter" for convenience) originates primarily in 
the optically thin outer disk, from several to perhaps hundreds of AU. The low column densities 
i n this region lead to a direct relationship between the submillimete r luminosity and the disk mass 



( Hildebrand 



1983 



Beckwith et al 



1990 



Andre &: Montmerld Il994l ) . Moreover, the shape of the 



long- wavelength SED is at l east partially controlled by the disk opacity (jBeckwith &: Sargent 



1991 



Mannings &: Emersonlll994l ). The optical properties of the d ust grain populatioi i j especially t he 
particle size distribution, are encoded in this latter parameter (jPollack et al.lll994l : lDrainell2006l ). 



In practice, one can infer the circumstellar dust mass around a young star with a submillimeter 
flux density and constrain the bulk optical characteristics of the solid particle population using 
submillimeter colors. Constraints on these properties are empirical prerequisites for understanding 
the mechanisms and timescales involved in the planet formation process. The progenitor disk 
mass not only sets a limit on the material composition and structure of a potential planetary 
system, but can also elucidate whether any particular mechanism of planet formation is more 
plausible than another. As an example, models that rely on gravitational instabilities in the disk 
to drive rapid giant pl anet formation are only tenable if the disk is sufficiently massive (iBoss 
19981 : iMayer et al.l 120071 ). Moreover, interpreting submillimeter colors in terms of the current solid 
particle size distribution in these disks may reveal clues on the efficiency of the grain growth process, 
representing the first tentative steps of building planetesimals. 



Beckwith et al.l (|l990l ) presented the first large survey of the 1.3 mm emission from young stellar 



objects (YSOs; defined here as a young star and its circumstellar en vironment) in the Taurus- A uriga 
star formation region. Following up with a suppl ementary survev (lOsterloh Beckwithlll995l) and 



smaller, multiwavelength submillimeter samples (jBeckwith &: Sargent 



1991 



Mannings &: Emerson 



19941 ). these studies provided the first detailed look at the circumstellar disk mass distribution 
and the shape of the opacity spectrum. The results argued for a typical disk mass similar to that 
required to form the Solar System and an opacity spectrum shallower t han noted in the inte r stellar 
medium (ISM), possibly due to grain growth. A round the same time, lAndre &: Montmerld (119941 ) 
built on their previous work (jAndre et al.lll99d ) with a pioneering investigation of the millimeter 
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continuum emission from YSOs in the p Ophiuchus dark clouds. One major impact of their work 
was highhghting that the milhmeter luminosities appeared tied to the same YSO evolution sequence 
typically tracked in the infrared: younger sources on this empirical sequence were shown to have 
larger circumstellar dust masses. 

To build on these landmark studies while taking advantage of technological advancements 
in the sensitivity and efficiency of submi llimeter detectors, we unde rtook a comprehensive, deep. 



multiwavelength survey of Taurus YSOs ([Andrews &: Williams! |2005| . hereafter Paper I). Centered 



around an 850 pm sample of ~150 sources, this work had the benefits of at least a factor of 5 
better sensitivity and a considerably more uniform luminosity limit than previous observations. In 
addition to a more complete census of disk masses and their evolution, this survey permitted us to 
derive the distribution of submillimeter colors between 350 /xm and 1.3 mm, and show that these 
also changed as a function of YSO evolutionary state. In this article, we present a similar study 
of the submillimeter continuum emission from circumstellar dust in the p Oph region. In §2 we 
describe the new survey data, its acquisition and reduction, and the compilation of an extended 
sample similar to the one used for Taurus. Disk masses and colors are computed in §3, and the 
results for p Oph YSOs are systematically compared with their counterparts in Taurus. Potential 
relationships of these properties with evolutionary diagnostics, stellar host parameters, and the 
planet formation process are discussed in §4, before concluding with a summary in §5. 



2. Observations and Data Reduction 



The Submillimeter Common-User Bolometer Array (SCUBA: [Holland et al.lll999l ) instrument 
at the 15 m James Clerk Maxwell Telescope (JCMT) was used to observe 48 p Oph YSOs during 
two runs in 2004 April and 2005 June. The observations were conducted in photometry mode, 
providing simultaneous measurements of the 450 and 850 /xm continuum emission within FWHM 
beam diameters of 9" and 14", respectively. Pointing centers accurate to < 1% of the beam size wer e 



determined from the Two Micron All Sky Survey (2MASS) Point Source Catalog (ICutri et al.ll2003l ). 
The typical precipitable water vapor level during the observations was < 1.8 mm, corresponding 
to zenith optical depths of 0.3—0.4 at 850 /im and ~2 at 450 pm.. The same observing and data 
reduction strategies described in detail in Paper I were utilized for this sample. The systematic 
uncertainties associated with absolute fiux calibration are expected to be ~10% and 25% for 850 
and 450 /im, respectively. 

Images of the 350 pm continuun i for 7 YSOs were o btained with the Submillimeter High Angu- 
lar Resolution Camera (SHARC-II: iDowell et al.l 120031 ) on the 10 m Caltech Submillimeter Obser- 
vatory (CSO) telescope in 2004 March and June. As with the SHARC-II observations described in 
Paper I, these diffraction-limited images (FWHM beam diameter of 9") were obtained by constantly 
scanning the telescope in the vicinity of the source in a Lissajous pattern. Further details on the 
observing technique and data reduction are highlighted in Paper I. The absolute fiux calibration of 
these 350 /im data is accurate to ~25%. 
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Unfortunately, the large submillimeter photometry survey of p Oph YSOs that was originally 
planned was cut short by the de-commissioning of SCUBA in the summer of 2006. For a more com- 
plete sample of p Oph similar in size and content to that presented in Paper I for Taurus, the data 
described abo ve were supplemented with a large number of 1.3 mm single-dish observations from 



1998 



Stanke et al 



the hterature (lAndre k, Montmerld Il994l : Ijensen et al.l 119961 : iNiirnberger et al.l 119981 : iMotte et al 



20061 ) ■ This expanded sample consists of 147 objects; 48 with multiwavelength 



data, and 99 others with only 1.3 mm flux densities from previous work. The vast majority of 
the sources lie in the L1688 dark cloud, with a few others scattered in L1689 and the surrounding 
Oph-Sco complex. Table [T] lists the 350, 450, 850, and 1300 p,m flux densities for these 147 sources, 
along with disk masses (§3.1), submillimeter colors (§3.2), and other relevant properties. Figured] 
summarizes some of the basic characteristics of the p Oph sample. 



In what follows, we follow IWilking et al.l (j2005l ) and adopt a distance of 150 pc to p Oph. 
For this distance, the projected FWHM beam sizes are 1350 AU for both 350 and 450 pm and 
2100 AU for 850 pva.. Because these single-dish beams are substantially larger than the circumstellar 



dust associated with Class II and III sources (e.g., iKitamura et al 



2002 



Andrews &: Williams 



20071 ). it is safe to assume that the values listed in Table [T] for such objects represent the total 



integrated continu um flux densities. As shown in mapping observations (e.g., iMotte et al 



1998 



Stanke et al.l l2006l ) and discussed in Paper I, this is not likely the case for Class I sources. It is 
important to emphasize that the photometry for the Class I YSOs listed in Table [1] corresponds 
to the combined continuum emission contributions from both a circumstellar disk and the inner 
regions of a circumstellar envelope. 



3. Results 



3.1. Disk Masses 



Following our work in Paper I, we aim to provide some empirical calibration of the relationship 
between submillimeter luminosities and disk masses through the use of a simple disk structure 
model. A geometrically thin irradiated dust disk emits thermal radiation modified by the optical 
properties of the solid particles from a continuous series of annuli at a range of stellocentric radii. 
The SED of this emission at a given frequency v (here defined as Ly = AttcPuF^) can be described 
in terms of the physical conditions in the disk as 

rRd 

Ly = STT^i^cosi / By{Tr) (1 - e-"'"-'-''"^') r dr, (1) 

where i is the disk inclination, tq and Ra the inner and outer radii, By{Tr) the Plan ck function at 
a radius-dependent temperature, and T,y_r the optical depth of the disk material (e.g., lAdams et al 



19871 : iBeckwith et al.lll990l : IChiang fc Goldreichlll997l l. The latter is taken to be the product of the 
opacity, k^-, and surface density, S^. The disk structure is approximated with radial power-laws, 
where oc and S.^ cx r~^. To simplify the SED modeling process, we define the opacity to have 
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a spatially homogeneous power-law spectrum such that oc with (3 = 1 and a normalization of 
kq = 0.1 cm^ at 1000 GHz (cf.. iBeckwith et al.lll99d ). This opacity implicitly assumes a 100:1 
mass ratio of gas to dust. 

We identified a subsample of 19 p Oph YSOs that had suitable mid-infrared through millimeter 
SED data in the literature to determine from a fit to this simple disk structure model. For 
consistency with the work done in Paper I for Taurus sources, we fixed the following parameters 
in the SED modeling of this subsample: i = 0°, rg = 0.01 AU, Rd = 100 AU, and p = 1.5. The 
SEDs longward of 5 pm. were de-reddened (see §4.2) and fit to the model using x^-minimization. 
The results of the SED fits are compiled in Table [21 including literature sources for the SED data, 
the reduced statistic, degrees of freedom in the fit, and the best-fit values for the 3 free model 
parameters; Ti (the temperature at 1 AU), q (the radial power-law index of the temperature profile), 
and Md (determined from the surface density profile and fixed outer radius). 

Using these results and their counterparts for the Taurus sample (see Paper I), Figure [2] 
demonstrates that the submillimeter fiux densities and best-fit disk masses are correlated. Most of 
the YSOs do not have sufficient data to determine from a SED fit, and therefore an alternative 
method is required. As dis cussed in the introdu ction, most disks are primarily optically thin at 
submillimeter wavelengths (IBeckwith et al.l Il990l ) . Assuming that r^^r ^ 1 and the emission is 
generated in a roughly isothermal region with a characteristic temperature (Tc), Eqn. 1 reduces to 



(2) 



This approximation explicitly shows that the submillimeter luminosity is a direct diagnostic of the 
total optical depth in the disk, or the product Kj^M^. The behavior of this simplified relationship 
is able to reproduce the SED fit results for both Taurus and p Oph YSOs quite well in Figure 2, 
where it is plotted as a solid line for Tc = 20 K. 

The disk masses (or upper limits) for the remaining 128 YSOs in this sample without sufficient 
SED data were computed according to Eqn. 2 with a characteristic temperature of 20 K. These 
Md values are compiled in Table [H with preference given to using the more modern and sensitive 
850 //m flux densities when available. The resulting cumulative distribution of p Oph disk masses for 
the full sample is shown in F igure [3l constructed with the Kaplan-Meier prod uct-limit estimator to 
incor porate the upper limits (IFeigelson &: Nelsonlll985l : see also e.g.. Paper I or I Andre &: Montmerk 
1994 ). This distribution (and others like it in what follows) shows the probability that a given disk 
in the sample has a mass greater than or equal to the abscissa value. The median disk mass in 
the fuh p Oph sample is {Md) « 0.002 Mq. Roughly 28% of these sources have Md > 0.01 Mq, a 
representative value for the minimal amount of solar-composition material require d to account for 
the p lanets in the Solar System (the Minimum Mass Solar Nebula, hereafter MMSN: IWeidenschilling 
19771 ). Only 40% of the p Oph YSOs have Md < 0.001 M©, corresponding to the mass of Jupiter. 
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3.2. Submillimeter Colors 



The submillimeter continuum spectra of YSOs generally follow a simple power-law behavior, 
such that Fi, oc iv". The spectral index a is essentially a submillimeter color, and potentially 
encapsulates important information about the optical properties of solid particles in the disk. 
For the optically thin, isothermal approximation (Eqn. 2), the submillimete r continuum einission 



has a spectral dependence Lp cx u^k^ oc in the Rayleigh-Jeans limit ( Beckwith et alJll99Cl : 



Beckwith &: Sargenti Il99ll ). With these assumptions, the submillimeter color is directly related 
to the shape of the opacity spectrum, a = 2 + /?, which in turn is thought to be a composite 
diagnostic of the dust particle size distribution, naineralqgical compos ition, and morphology (e.g., 
Mivake &: Nakagawal Il993l : iHenning &: Stognienkd 119961 : lDraindl2006l ). In the ISM, the generally 
small (sub-rnicroii ) dust particles are relatively inefficient submillimeter emitters with P ~ 1.7 (e.g .. 



Li &: Draind l200lh . resulting in steep submillimeter spectra with a ~ 4 (e.g., iHildebrandl Il983l ). 
However, in the high-density environment of a circumstellar disk, solid particles are e xpected to grow 



to su bstantially larger sizes via collisional agglomeration (see the recent review by iDominik et al 



20071 ). For a given grain composi tion and morphology, increased particle sizes due to such growth 



will decrease the value of (3 (e.g., iMiyake &: Nakagawa 



1993 



Draind 120061 ). This in turn would be 



manifested as a shallow submillimeter continuum spectrum, or decreased a. 

To assess the potential signatures of grain growth in p Oph disks, we have calculated submil- 
limeter colors (a) between 350 /^m and 1.3 mm whenever multiwavelength data were available. In 
cases with >2 flux density measurements in this wavelength range, a was determined by a linear 
fit in the {logz^, logi^i^} plane. Otherwise, a simple flux ratio was utilized to measure a color or 
upper limit. The inferred values of a are included in Tabled! Figure [4] shows their cumulative dis- 
tribution, again constructed with the Kaplan-Meier estimator to include upper limits. The median 
value is (a) ~ 2, with only a few percent of the YSOs showing a > 3. 

Based on the above discussion, it would be tempting to conclude that the submillimeter SEDs 
for the p Oph sample show clear evidence for significant solid particle growth (with /3 ^ 1). However, 
there are a number of other effects that can serve to mimic the signatures of a shallow opacity 
spectrum. Perhaps the most relevant of these is optically thick contamination, originating in the 
dense inner disk. As a tracer of the disk temperature, emission with large optical depths will bias 
a to lower values because it has the same spectral shape as the Planck function in the Rayleigh- 
Jeans limit (a = 2). Interferometric observations that resolve the disk structure can be used to 
quantify the c ontamination, and usually indicate optically thick f ractions of < 20 — 30% shortward 



of 1 mm (e.g., iTesti et al 



2001 



2003 



Andrews &: Williama 120071 ). Moreover, the presumed origin 



of the submillimeter continu um emission near the cold disk midplane (e.g., IChiang &: Goldreich 
19971 : iDullemond et al.ll2002l ) implies that the Rayleigh-Jeans limit may not be valid (i.e., if much 
of the emission comes from dust with T < 15 K). In this case, the thermal emission is considerably 
shallower than and may contribute to the observed low a values. Finally, there are more technical 
considerations, particularly the rather large uncertainties on a due to calibration inaccuracies and 
small frequency leverage, as well as contamination from non-disk emission (see Paper I for details). 
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Despite these uncertainties, it would be difficult to account for the sub millimeter color distr i- 
bution in Figure H] for an opacity spectrum similar to that in the ISM ^see b'Alessio et alJ[200ll ). 
Detailed disk structure modeling and spatially resolved data are required to reach a more definitive 
conclusion. In those cases where the disk structure has been resolved and can be u sed to correct the 
color s (a) for optically thic k contaminatioii , the t ypical opacity index is /3 ~ 1 (|Testi et al.l 12001 



2003; 



Rodmann et al 



2006; 



Lommen et al. 



2007; 



Andrews Williama 120071 ) . Without resolved 



data, the colors presented here serve solely as empirical properties of YSOs. While the opacity 
spectra of these sources remain uncertain, the observed colors are certainly significantly different 
than noted for dust in the ISM, in the same sense as expected for when particle sizes are increased. 
The issue of grain growth will be addressed further in §4.3. 



3.3. A Comparison with Disks in Taurus- Auriga 



Having established the basic submillimeter properties of YSOs in p Oph, we can make a com- 
parison with their counterparts in Taurus from Paper I. Taurus is a s ubstantially less cr owded 
region, with projected stellar densities reaching only ^50—6 stars pc~^ (Goinez et al.lll993l ) com- 
pared with the ~200 stars pc~^ noted in the p Oph core ( Greene &: Young] 1992 ). Despite the 
higher concentration of YSOs in p Oph, a simple calculation of the collision timescale shows that 
dyn amical interactions are not expected to significantly affect their disk properties. De-projecting 
the iGreene Young] ( 19921 ) stellar density into a 1 pc^ volume (n* ~ 200 pc~^) and assuming a 



typical disk size (r ~ 100 AU =^ a ~ 10 ^ pc ^, where a is the coUisional cross section) and cluster 



velocity dispersion ~ 1 km s ^) implies a disk— disk collision timescale of tc 
roughly 3 orders of magnitude larger than the cluster lifetime. 



IGyr, 



In addition to being more spatially concentrated, p Oph members are often still deeply em- 
bedded in their natal molecular cloud material. This is in stark contrast to most YSOs in the 
Taurus association, which by comparison are relatively unobscured. Given these environmental 
distinctions, it is worthwhile to examine if they have any impact on disk properties by directly 
comparing the results for the two regions. 

To make the comparison a fair one, most of the Taurus disk masses were re-calculated from 
the 850 pm flux densities in Paper I using Eqn. 2, the same parameters adopted here (see §3.1), and 
the appropriate distance, d = 140 pc. A subsample of Taurus sources have sufficient information to 
determine from a fit to the full SED (see Table 2 in Paper I), and in those cases the best-fit disk 
masses were adopted. Direct comparisons of the cumulative distribution functions of and a for 
the similar-sized (~150 YSOs) p Oph and Taurus samples are shown in the top panels of Figure 
O These distribution functions take into account upper limits using the Kaplan-Meier estimator 
and should be a reliable comparison for > 0.001 Mq, corresponding to a rough luminosity limit 
for the p Oph sample. The error envelopes of the Taurus distributions are shown in grayscale and 
the p Oph distributions as histogram s with error bars. The tw o-sample statistical tests described 
in the survival analysis formalism of iFeigelson NelsonI (Il985l . i.e., the Gehan, Logrank, Peto & 
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Peto, and Peto & Prentice tests) confirm what is apparent from Figure [5j the distributions of 
and a for p Oph and Taurus YSOs are essentially indistinguishable. The results of these tests are 
provided in the top half of Table [3l 

Since the relative numbers of YSOs in various evolutionary states are not the same in the p 
Oph and Taurus samples, a closer look at more similar subsamples makes a better comparison. 
The shap e of the infrared SEP is often utilized as an empirical diagn ostic of the YSO evolution 
sequence dLada fc Wilkinglll984l : kdams k Slu3ll98d : kdams et al1ll987l l . The parameter of interest 
is the SED ir idex n, defined froin ~2— 25 pm via i^Fi, oc v"". We adopt the classification breakdown 
suggested bv lGreene et al.l (jl994l ). where Class I objects (star+disk+envelope) have n < —0.3, Flat- 
Spectrum objects (star+disk+envelope?) have —0.3 < n < 0.3, Class II objects (star+disk) have 
0.3 < n < 1.6, and Class III objects (star only) have n > 1.6 . Infrared photornetry for the p Oph 



samp l e was collected from the literature (l\yilking et al 



1994 



Bontemps et al 



2001 



Barsony et al 



2005 



198S: 



Weaver Jones 



1992 



Greene et al 



Padgett et al.l l2006l ;i. the 2MASS Point Source 



Catalog (jCutri et al.M2003l ). and the preliminary results of the "Cores to Disks" Spitzer Legacy 
Survey of the region (L. Allen 2007, private communication). The infrared photometry of Taurus 
sources utilized f or this sam e purpose in Pa p er I w as supplemented with Spitzer photometry from 
Hartmann et al.l (j2005l ) and iLuhman et al.l (j2006l ). Infrared SED slopes (n) were computed for 
each YSO in the p Oph and Taurus samples with a simple power-law fit to the SED, corrected 
for reddening whenever possible (extinctions were computed as described in §4.2 or Paper I). The 
resulting SED classifications for the p Oph sample are included in Table [TJ 

This empirical evolution sequence was used to generate YSO subsamples in both star-forming 
environments and compare them with the aforementioned two-sample statistical tests. The results 
are compiled in the top portion of Table O and the cumulative distributions of and a for 
various SED classifications are also shown in Figure [5j The P values in Table [3] correspond to the 
probability that the two subsamples are dra wn from different parent distributions. The ranges of 
P represent the various statistical tests (see iFeigelson &: NelsonI Il985l ) . After concluding that the 
few Flat-Spectrum sources do not substantially influence the results (as was noted in Paper I), they 
were incorporated into the Class II subsample. 

Although not significant by formal statistical arguments, the distributions of and a for 
Class I objects in the two regions appear somewhat different in Figure [5j Little import should be 
attached to this apparent discrepancy, as the photometry has not been conducted in a systematic 
way and the potential for different amounts of envelope contamination has therefore not been 
properly taken into account. On the other hand, the level of detail in the agreement between 
the p Oph and Taurus distributions of and a for Class II sources is remarkable. For disk 
masses greater than ~0.001 M©, the distributions between the two regions are essentially identical, 
even showing the same dip feature around the MMSN value (0.01 M©) and decrease to virtually 
zero probability for disks with > 0.1 Mq. The median Class II disk in both samples has 
(Md) ~ 0.005 Mq, with roughly one third of each sample harboring disks with > Mmmsn- 
The shapes of the submillimeter color distributions between the two regions are also fairly precise 
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matches, with median values {a) ~ 2 and only a few percent of each sample having a > 3. 

As with the full samples, the subsamples described above indicate that there are no statistically 
significant differences between the submillimeter properties of YSOs in the p Oph and Taurus 
regions at any given evolutionary state. Given the roughly equivalent sizes and basic characteristics 
of the samples in these two regions and the homogeneous methodology utilized to infer and a, 
this result clearly indicates that the local environment of a low-mass, low stellar density (<C 10^ stars 
pc~^) region does not substantially affect the submillimeter properties of the associated disks. 



4. Discussion 



4.1. Disk Evolution 



While confirming the similarities between the submillimeter properties of circumstellar dust in 
the p Oph and Taurus regions, the panels in Figure [5] also show that there are substantial changes 
in Md and a across the empirical evolution sequence characterized by the infrared SED index. 
Two-sample tests on the censored datasets directly confirm that the distributions of Class I, II, 
and III sources in p Oph are statistically different. The test results are catalogued in the bottom 
half of Table [3l As with the YSOs in the Taurus clouds (Paper I), submillimeter luminosities (and 
therefore disk masses) and colors decrease substantially along the Class I — > II — > III sequence. 

Figure [6] directly demonstrates how the disk masses and submillimeter colors vary as a function 
of the infrared SED index (n), now combining the data for the p Oph and Taurus samples. The 
top panel reinforces what is seen in the left panels of Figure Md decreases as the infrared 
SED becomes bluer (i.e., n increases) along the evolution sequence. Using only the detections, a 
correlation between and n is present at the 4 a level, with a Spearman rank correlation coefficient 
of —0.32. The correlation significance is enhan ced when th e uppe r limits are incorporated with any 
of the survival analysis statistics discussed bv llsobe et al.l (|l986l ). Disk masses decrease smoothly 
with a shallow slope for YSOs with infrared SEDs redward of the Class II/III boundary (at n ~ 1.6), 
but then drop precipitously. This no n-linear behavior in the {n, M^j-plane was in fact originally 
postulated in the pioneering study bv I Andre &: Montmerld (| 19941 ) . along with the first indications 
that the p Oph and Taurus disk mass distributions were the same. In an effort to approximately 
determine the depth of the drop-off beyond the Class III boundary, we stacked the 850 pm SCUBA 
photometry observations of all undetected Cla ss III sources in p Oph and Taurus (Paper I) to 
derive their average luminosity value (see, e.g., ICarpented l2002l ) . The mean 850 /xm fiux density 
for the 54 stacked Class III non-detections (5 in p Oph and 49 in Taurus) was 0.1 it 0.3 mJy. With 
the same assumptions adopted in §3.1, the 3 a "average" limit (0.9 mJy) indicates that the typical 
undetected Class III source has < 5 x 10^^ M©, roughly the mass of Neptune. 

The bottom panel in Figure [6] examines the comparative evolution of colors in the infrared 
and submillimeter bands. There is a marginal (2.9 a) correlation between a and n in this diagram 



- 10 - 



with a Spearman rank correlation coefficient of —0.30, suggesting tliat a systematically decreases 
along the infrared evolution sequence. A similar correlation for the Taurus sources alone was noted 
in Paper I, and there partially attributed to a genuinely evolving disk opacity. The idea behind 
this argument is that disks are expected to become more optically thin as YSOs evolve along the 
infrared SED sequence. In the case of no opacity evolution, a decreasing fraction of optically 
thick submillimeter emission as n increases should lead to increasing a, the opposite of what is 
noted in Figure [H On the other hand, the growth of solid particles as the circumstellar dust 
evolves along the sequence would result in a flattening of the opacity spectrum (decreased and 
therefore the observed decreasing a. An alternative (or possibly cooperative) scenario could be a 
decreasing tem perature in the disk m idplane (e.g., due to gravitational sedimentation of dust to the 
disk midplane: Id 'Alessio et al.ll2006l ). These possibilities can potentially be explored with spatially 
resolved multiwavelength data, although with such a weak correlation many more sources would 
need to be observed with interferometers. 

The above analysis confirms for p Oph what was shown for Taurus in Paper I: the infrared and 
submillimeter properties of circumstellar dust around young stars evolve along a similar empirical 
sequence. Given where such emission originates, this result suggests that the material in the inner 
(infrared) and outer (submillimeter) parts of disks is significantly affected by global changes in disk 
properties. One of the obvious candidates for the cause of such evolution is the metamorphosis 
of circumstellar structure, including the dissipation of the envelope (Class I — > II), the viscous 
evolution of the accretion disk (through Class II), and the dissipation of the disk (Class II — > III). 
A perhaps equally important contribution to this evolution could be the result of the substantial 
growth of solid particles. In reality, the submillimeter luminosity is determined by the product 
KyM^: even if the total circumstellar mass remained the same across the evolution sequence, the 
observed trends could both be explained by the changes in the opacity expected from grain growth. 

Given the abruptness of the transition at the Class II/III boundary in the top panel of Figure [6] 
and the low submillimeter detection rate for Class III sources (see also Paper I) , it appears that the 
signatures of the inner and outer disk (infrared and submillimeter excesses, respectively) disappear 
within a short time of one another. In Paper I, we used the standard "duty-cycle" argument for this 
based on the very small percentage (< 3%) of YSOs that have long- wavelength (i.e., submillimeter) 
emission from the outer disk, but no near-infrared excess from the inner disk. The scarcity of 
such "transition" objects suggests that if all disks go through such a phase, it must be rather 
short in duration (< 10^ yr). A vari ety of mechanisms have been proposed for this rapid radial 
evolution, from planet form ation fseelNaiita et al.ll2007l. and referenc es therein) to rapid mass loss 



in photoevaporative winds (j Clarke et al 



2001 



Alexander et al.l 12009 ). 



4.2. Relation to Stellar Properties 



Despite the above effort to associate changes in the submillimeter properties of circumstellar 
dust with an empirical evolution sequence, it is still not possible to tie this sequence to an abso- 
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lute timeframe. An alternative examination of disk evolution can be made by relating disk and 
stellar parameters, particularly masses (M*) and ages (t*). To determine the latter, we first col- 
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their stellar masses and ages. For the remainder of the sample, we used the s pectral type— effective 
tempe rature scale, intrinsic colors, and bolometric corrections tabulated bv lKenvon Hartmann 
(jl995l ). Visua l extinctions were c alculated from color excesses (in the optical-red whenever avail- 
able) and the ICohen et al.l (|l98ll ) extinction law. Stellar luminosities were then computed from 
the de-reddened photometry a nd appropriate bolomet r ic cor rection. Stellar masses and ages were 
determined by reference to the lD'Antona &: Mazzitellil (|l997l ) theoretical pre— main-sequence mass 
tracks and isochrones in an H-R diagram. 

This is the sam e technique ad o pted f or the Taurus YSOs in Paper I, and is similar in practice 



to that utilized by IWilking et al.l (|2005l ). The stellar mass and age distributions were already 
highlighted in Figure [TJ The p Oph sample is slightly younger than the Taurus sample, with a 
median age of ~0.7Myr compared to ~lMyr. The cumulative distributions of disk-to-star mass 
ratios (M^/M*) for Class II sources in both the p Oph (lines) and Taurus (grayscale) samples are 
shown together in Figure [71 The p Oph sources generally exhibit slightly higher mass ratios than 
their Taurus counterparts, although the M^/M^ distributions for the two regions are not formally 
different in a statistical sense. The median ratio in p Oph is {Md/M^) ^ 1.0%, compared to ~0.8% 
for Taurus. In this case, selection effects may be mimicking a subtle diagnostic of disk evolution 
from the slightly younger p Oph sample to the more mature Taurus YSOs. The p Oph sample does 
have a somewhat different distribution, suggesting that any true hint of a change in M^/M^, as 
a function of cluster age should await larger samples with more similarity in their stellar properties. 

Figure El exhibits the variations of M^, M^/M,,, and a with time for the combined p Oph and 
Taurus sample. While there is a large scatter of disk properties at any given time, some marginal 
trends are apparent. Although a definitive correlation between and remains unconfirmed, 
there are clearly more non-detections as the stars get older. Using survival analysis techniques 
to incorporate upper limits, the Cox proportional hazard model, generalized Kendall's tau, and 
Spearman's rho correlation tests (|lsobe et al.l 119861 ) all suggest that the disk-to-star mass ratio is 
indeed correlated with the stellar age, dropping from a few percent at ~0.1Myr to a factor of 30 
or more lower by 10 Myr. It would be useful to measure disk masses for clusters with slightly older 
ages to confirm such trends, so that the strength of the correlation does not hinge on the youngest 
sources (i.e., t < 10^ yr), where the isochrones are less trustworthy. No obvious trend between a 
and age is apparent. 



Before moving on, it is worthwhile to comment on the adopted distance {d = 150 pc) and how 
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alternative values would be expected to affect the above analysis. Following IWilking et alj (|2005l ) 
and others, we based the choice of d on the typical range of '^130— 170 pc from m easurements of 
stars in the larger Upper Scorpius complex that encompasses th e p Oph clouds ( e.g., IWhittetlll974l : 



de Geus et alJ Il989l : iKnude Sz Hod Il998l : Ide Zeeuw et alJ Il999l : iMamaiekl 120071 ') . However, other 
authors have focused on the near end of that range, using the d = 125 pc advocated by de Geus and 
colleagues based on two arguments suggesting that the dark clouds are in the Upper Sco foreground: 
(a) the morpholog y of the distance mo dulus vs. extinction diagram for Upper Sco sources near the 
p Oph cloud edge (Ide Geus et al.lll989l). an d (6) some circumstantial evidence from the morphology 
of the local ISM (jde Geus k Burtonlll992l l. 



If this closer distance is adopted, all luminosities utilized above are reduced by ~30%, with two 
important consequences. First, the decreased bolometric luminosities in the H-R diagram would 
result in considerably older ages for the p Oph sample (stellar masses would remain approximately 
the same, owing to the essentially vertical tracks in this diagram). With the roughly logarithmic 
behavior of isochrones along the H-R diagram luminosity axis, the increase in could be a factor of 
2—4. Second, this nearer distance would lead to a substantial discrepancy in the p Oph and Taurus 
Mfi distributions, such that disks in the former region would be significantly less massive. Taken 
together, using a distance as close as 125 pc would imply that p Oph YSOs are older than their 
Taurus counterparts, and that disk masses have decreased substantially over the corresponding age 
difference. Obviously using ages from the H-R diagram is not an ideal way to constrain a cluster 
distance. However, in this case it does not really make sense that the much more deeply embedded 
YSOs in p Oph are older than their essentially unobscured Taurus counterparts. Perhaps this issue 
can be resolved with a more complete census of th e ste llar population in the p Oph clouds, along 
the lines of the work done by IWilking et al.l (j2005l ) and iDoppmann et al.l ([2003, l2005l ) . 



4.3. Benchmark Submillimeter Properties of Circumstellar Disks 

Now that large submillimeter photometry samples of YSOs for the p Oph and Taurus associa- 
tions have been collected, it is appropriate to present quantitative reference values for some of the 
basic properties of their circumstellar environments. The cumulative distributions of the censored 
datasets presented above were used to determine the median values of M^, M^/M,,, and a com- 
piled in Table |4] for various evolutionary states. Here we have distinguished the Class HI sources 
(n > 1.6) into two categories: those with submillimeter detections are termed "transition" objects 
(because dust in the inner disk is largely absent, while a significant outer disk remains) and those 
undetected sources retain the standard Class HI label. These benchmark values clearly highlight 
the evolution in the submillimeter properties that have been carefully determined in this work and 
Paper I. Figure [9] graphically illustrates the trend of decreasing (M^) along this sequence, with 
representative error bars showing the distribution quartiles. Typical disk masses are similar to the 
MMSN for much of their evolution, only dropping substantially at the end of the Class II stage. 
The transition disks have a typical disk mass roughly an order of magnitude lower than the MMSN 
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value, corresponding roughly to a Jupiter mass of material. The stringent stacked upper limit 
placed on the remaining Class III disks is more than an order of magnitude lower still, suggesting 
that primordial circumstellar disks at this stage in evolution are probably non-existent. However, 
this limit still does not r ule out the presence of low-mass debris disks similar to those noted around 
older nearby stars (e.g.. [Holland et al.lll998l ). 



The accretion disks present in the Class H stage of evolution are generally regarded as the 
birthsites of planetary systems. As such, their properties are directly relevant for constraining 
the initial conditions available for the planet formation process. In light of their significance in 
this matter. Figure [10] shows the differential and a distributions for the p Oph and Taurus 
Class II sources that have firm submillimeter detections. Compiled from 125 individual sources, 
this distribution of disk masses peaks at ~0.01 Mq (the MMSN value) and covers a fairly wide 
dispersion. As would be expected from the analysis of individual regions presented here and in 
Paper I, the a distribution for these sources has a relatively narrow peak centered around a ~ 2. 
With the similarities between the disk properties in the two star-forming regions (§3.3), these values 
should represent a combined, universal description of the most basic submillimeter properties of 
Class II disks. 

But how precisely are the values of these properties constrained in an absolute sense? Of 
course, the submillimeter colors are empirically measured quantities and, aside from occasionally 
large error bars, their values are firmly determined. The challenge is interpreting the a values 
in terms of the solid particle population in these disks. The evidence is now fairly clear on the 
observational side that the typical disk opacity spectrum is more shallow than in the ISM, with Bjy 1 



deemed appropriate to explain the submillimeter SEDs (e.g.. Paper I: iBeckwith k. Sargent 



1991 



Mannings &: Emersonlll994l : lRodmann et al.ll2006l : [Andrews &: Williamsll2007l ) . This effect has been 
interpreted as the signature of an increase in the maxim um dust particle size, frorn ~sub-micron 



sizes in the ISM to ^millirneter sizes in a typical disk (iMiva ke & Nakagawa 



1994 



Henning Stognienko 



1996 



D'Alessio et al. 



2001 



1993 



Pollack et al 



Drain&.2006i ). This amounts to ~4 orders 



of magnitude in particle growth from the time when these dust gra i ns we re incorporated into 
circumstellar disks/envelopes from the ISM. Moreover. iD'Alessio et al.l (120011 ) emphasize that this 
actually represents a lower limit on the amount of particle growth, as (3 remains near unity for 
maximum particle sizes well beyond the (sub-)millimeter wavelengths used to constrain it (see 
their Fig. 3). In essence, submillimeter colors indicate that significant grain growth has occurred 
in these disks, although it remains unclear exactly how far along it may have proceeded beyond 
millimeter sizes. 

In a sense, the opacity acts as an accounting system for the mass budget in a disk. Not all solid 
particle sizes are generating sufficiently detectable emission at submillimeter wavelengths, and the 
opacity must also represent the mass contribution of the inefficient emitters in the conversion of a 
submillimeter luminosity to a disk mass. A direct, non-degenerate measurement of in a disk is 
not currently possible, and so a standard approximat e value has been fixed and regularly utilized 
in this field. In their admittedly simplified analysis, ID'Alessio et al.l (j200ll ) show that this fixed 
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opacity value, Ki.smm ~ 0.02 cm^ (as well as /3 = 1), is appropriate when the largest particle 
sizes are ~1 mm. However, the opacity is significantly diminished when the maximum particle size 
is larger than this value. Because oc Lsmm/i^v, an underestimate of the amount of particle 
growth in these disks implies an underestimate of their masses. 

Using labora tory experimen t s and computer models as a guide for the collision conditions in 
the disk interior, iDominik et al.l (|2007l ) argue tha t solid bodies u p to ^^ meter size scales can be 
formed on rather short timescales. The work by iD'Alessio et al.l (|200ll ) shows that meter-sized 
bodies would have an opacity ~10x lower than what has been assumed here (but with the same 
value of /3) . If such growth is typical, the inferred disk masses may be underestimated by an order 
of magnitude. iHartmann et al.l (|l998l ) proposed an interesting test on the validity of the standard 
opacity prescription. Aside from the method used here based on the submillimeter continuum 
emission from dust, they argued that mass accretion rates (M) determined from optical/ultraviolet 
spectra can be used to estimate with some reasonable assumptions about the accretion history. 
Their disk mass estimates, « 2t*M, are based solely on observations of the gas phase in the 



disks , and are therefore completely independent of the values presented here (see IHartmann et al 



19981 . regarding the numerical factor of 2). Using the stella r age s derived in t he pr evious section 
and the accretion rates compiled by IHartmann et al.l (|l998l ) and iNatta et al.l (|2006l ) , disk masses 
were calculated using both methods for the subset of p Oph and Taurus Class II sources with the 
required information available (i.e., M, t*, and a millimeter detection) 

A comparison of the two methods for computing in Figure [TT] clearly demonstrates that 
disk masses estimated from dust emission are systematically smaller than those inferred from the 
accretion rates. The dotted line in this plot has a slope of unity and bisects the sample into equal- 
sized groups: it lies roughly an order of magnitude below the solid line marking equality between 
the two Mrf determination methods. Of course, both methods of estimating involve significant 
uncertainties. The unknown disk structure in most cases can lead to relatively small uncertainties 
on the values determined from submillimeter continuum emission (s ee Paper I; And r ews fc 
Williams 2007). Larger uncertainties plague the estimates from the Hartmann et al. (jl998l ) 
method, due to both the subtleties of determining M and the challenges of specifying individual 
stellar ages. While these undoubtedly contribute to the scatter in Figure \TT\ they are unlikely to 
cause the noted systematic offset. For example, a temperature uncertainty could easily translate 
to a 50% Md uncertainty (see Eqn. 2) for a given disk. But, the scatter in Figure [2] shows that it 
is improbable that the average adopted Tc is systematically 2x lower than assumed (let alone the 
10 X needed to explain Figure [TT]) . A similar argument holds for other disk structure parameters 
and the stellar ages used in determining from the accretion rate. 

Perhaps a variety of uncertainties can unfortunately conspire to produce the observed offset. 
To be fair, this comparison of estimates from different material phases and spatial locations 



^We preferred to re-calculate the accretion rates provided by Natta et al. ( 20od ) based on the stellar parameters 
adopted throughout this paper. Regardless, the results in Figure [TT] stand irrespective of whose numbers are used. 
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(gas near the inner rim and dust spread to the outer radius) rehes on the assumption that accretion 
disks can be approximated with relatively simple physical prescriptions. However, given all of these 
caveats, Figure [TT] is really the only check on the absolute precision of measurements currently 
available. If solely physical explanations are considered to explain these results, the growth of solid 
particles to centimeter or meter size scales (rather than the millimeter scales that correspond to 
the standard opacity prescription) is perhaps the most likely cause. Such growth is observationally 
reinforced by the rem arkable centimeter- wave dust emission from large particles recently detected 
for the TW Hya disk (|wilner et al.lbonfil ). 



As discussed above, the decrease in the opacity expected from this level of grain growth cor- 
responds to disk masses being under-estimated by up to an order of magnitude. If this is the case, 
the implications for the planet formation process are significant. The appropriate shift in the Md 
distributions presented here to account for the over-estimate could imply that a relatively large 
fraction (perhaps up to one third) of ^ ^1 Myr old Clas s II disks are massive enough to be gravita- 
tionally unstable (i.e., M^/M* > 0.2; IShu et al.lll99Cll ). If these instabilities are indeed common, 
the giant planet formatioi i process inay be rapidly accel erated via direct condensations from spiral 



density waves/rings (e.g.. Boss 
and gas capture models (e.g., 



1998; 



Maver et al 



Pollack et al 



1996 



20071) ■ compared to the standard core accretion 



Hubickvi et al.ll2005l ). A concerted effort utiliz- 



ing both high spatial and spectral resolution multiwavelength submillimeter data (e.g., from the 
ALMA and e-VLA interferometers) and state-of-the-art disk models may help shed some light on 
this subject in the near future. 



5. Summary 

We have presented new multiwavelength submillimeter continuum observations of 48 YSOs in 
the p Oph dark clouds. Supplementing this survey with 99 millimeter measurements in the litera- 
ture, we study the circumstellar dust around young stars in this region; the results are summarized 
here. 

• The statistical distributions of disk masses (Af^), disk-to-star mass ratios, and submillimeter 
colors (a; where Fy (x u") were derive d and directly compared with those for the Taurus 



association (j Andrews fc Williamsl 120051 ). Within the uncertainties, there are no significant 



differences between the submillimeter properties of YSOs in the two regions. 

There are statistically significant decreases in and a along the canonical YSO evolution 
sequence based on the shape of the infrared SED. Changes in the structure of the circumstellar 
environment and/or the average particle properties may explain the correlations. 

A tentative correlation between the disk-to-star mass ratio and stellar age are noted for the 
combined samples in both regions. Confirmation with older sources would be beneficial to 
assess whether or not this relationship depends on isochrone models at early ages (~ 10^ yr). 
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The Class II YSOs in p Oph and Taurus are utilized to establish benchmark values for sub- 
millimeter properties in the population considered to be the birthsites of planetary systems: 
Md ~ 0.005 M0, Md/M^ ~ 1%, and a ~ 2. 



An independent disk mass estimate based on the accretion rate (jHartmann et al.lll998l ) is used 
to compare with the submillimeter continuum measurements. The latter are systematically 
smaller than the former by roughly an order of magnitude on average. This may be the result 
of substantial particle growth (up to ~meter sizes), leading to an over-estimate of the opacity 
used to compute M^, from the submillimeter luminosity. Such disk mass underestimates may 
imply that some fraction of Class II disks are marginally gravitationally unstable, which could 
have important implications for disk evolution and giant planet formation. 
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Table 1. Submillimeter Properties of Circumstellar Disks in p Ophiuchus 













TP /O Cn I ^■rv-.^ 


ri^y\..6 mmj 








(JDject 




op i 


[mjy] 


[mJy] 


[mJy] 


[mjyj 




a 


"\T ^4- j^^-i 

IN otes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(8) 


(9) 


(10) 


Wa Oph 1, HBC 630 


III 


K2 




< 114 


< 8 


< 20 


< 0.0005 




1 


AS 205, HBC 254 


n 


K5 


4683 ± 118 


3280 ± 271 


891 ± 11 


450 ± 10 


0.03 


1.85 ± 0.21 


1 


Wa Oph 2, HBC 633 


HI 


Kl 




< 129 


< 9 


< 20 


< 0.0006 




1 


Wa Oph 3, HBC 634 


III 


KO 




< 145 


< 10 


< 25 


< 0.0006 




1 


WSB 3 


III 


M3 






< 34 




< 0.002 




2 


WSB 4 


HI 


M3 








< 13 


< 0.003 




3 


WSB 11 


H 








< 26 




< 0.001 




2 


L1719 B, 16191-1936 


FS 










220 ± 10 


0.04 




1 


WSB 18 


HI 


M2 






< 46 




< 0.003 




2 


WSB 19 


H 


M3 






< 78 




< 0.004 




2 


DoAr 16, HBC 257 


H 


K6 




< 361 


47 ± 8 


< 50 


0.003 


< 3.21 


1 


SR 1, EL 9, GSS 9 


HI 


B2 




< 154 


< 11 


< 50 


< 0.0007 




1 


SR 22, DoAr 19 


II 


M4 




< 162 


31 ± 3 


< 20 


0.002 


< 2.60 


1 


ROXs 2, VSS 20 


HI 


MO 




< 156 


< 11 


< 30 


< 0.0007 




1 


IRS 2, YLW 19 


H 


K3 








< 25 


< 0.005 




1 


IRS 7 


II 










< 40 


< 0.008 




1 


IRS 8 


HI 










< 40 


< 0.008 




1 


IRS 9, YLW 23 


HI 










< 10 


< 0.002 




1 


ROXs 3, HBC 636 


HI 


MO 








< 30 


< 0.006 




1 


IRS 10, ROXs 4 


HI 


K6 








< 28 


< 0.005 




3 


ROXs 5 


HI 


K7 








< 19 


< 0.004 




3 


SR 4, AS 206 


H 


K5 




< 768 


142 ± 7 


31 ± 6 


0.004 


1.50 ± 0.62 


3 


IRS 13, ROXs 7 


HI 


K6 








< 20 


< 0.004 




1 


DoAr 21, HBC 637 


HI 


Kl 








< 30 


< 0.006 




4 


SR 3, IRS 16 


HI 


AO 








< 25 


< 0.005 




1 


GSS 26 


II 


K8 




650 ± 123 


298 ± 7 


125 ± 20 


0.04 


1.65 ± 0.36 


4 


EL 18, GSS 29 


II 


K6 








< 10 


< 0.002 




4 


DoAr 24 HBC 638 


II 


K5 








< 30 


< 0.006 




4 


CRBR 12, ISO-Oph 21 


I 










85 ±5 


0.02 




4 


DoAr 23, WSB 26 


HI 


MO 






< 77 




< 0.004 




2 


EL 20, VSSG 1 


II 






< 666 


240 ±7 


50 ±10 


0.03 


3.69 ± 0.71 


1 


CRBR 15 


I 


M5 








40 ±5 


0.008 




4 


EL 21, GSS 30 


I 






1403 ± 116 


222 ± 5 


90 ±30 


0.01 


2.50 ±0.33 


4 


LFAM 1 


FS 






2662 ± 219 


572 ±8 


250 ± 20 


0.04 


2.18 ±0.31 


4 


GY 10 


HI 


M9 








< 10 


< 0.002 




4 


GY 11 


FS 


M7 








< 10 


< 0.002 




4 


DoAr 24E, GSS 31 


II 


Kl 




779 ± 121 


158 ±6 


70 ±20 


0.008 


2.32 ± 0.43 


4 


DoAr 25, WSB 29 


II 


K5 






461 ± 11 


280 ± 10 


0.03 


1.17 ±0.54 


1 


LFAM 3, GY 21 


FS 










110 ± 10 


0.02 




4 


ROXs 9A 


HI 


MO 








< 22 


< 0.004 




3 


EL 24, WSB 31 


II 


K6 


5242 ± 213 


3012 ±216 


838 ±8 


390 ± 10 


0.2 


1.99 ±0.22 


1 


ROXs 12 


HI 


MO 








< 19 


< 0.004 




3 


VSSG 27, GY 51 


II 










< 20 


< 0.004 




4 


ROXs 9B 


HI 


K2 








< 11 


< 0.002 




3 
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Table 1 — Continued 















r .o 111111 1 








V_y U J CL-lj 


SED 


SpT 


fm Tvl 


\m Tvl 


[mjyj 


fm Tvl 
[mjyj 


1™©1 






\'^) 












y ) 


\°) 






IRS 14, GY 54 


II 










< 30 


< 0.006 




1 


IRS 15, GY 58 


II 










< 30 


< 0.006 




1 


SI, EL 25, GSS 35 


HI 


B3 








< 10 


< 0.002 




4 


ROXs 9C 


III 


K4 








< 14 


< 0.003 




3 


GY 91 


I 










300 ± 15 


0.06 




4 


WSB 37, GY 93 


II 


M5 








< 25 


< 0.005 




1 


WL 8, GY 96 


II 










< 10 


< 0.002 




4 


IRS 17, WL 7 


FS 










< 30 


< 0.006 




4 


ROC 16, GY 101 


FS 










< 15 


< 0.003 




4 


EL 26, GSS 37 


II 


MO 






< 51 


15 ±5 


0.005 


< 2.88 


1 


IRS 19, VSSG 11 


III 


MO 








< 15 


< 0.003 




1 


WL 12, YLW 2 


I 






749 ± 130 


216 ±7 


115 ± 10 


0.01 


2.00 ± 0.36 


4 


EL 27, GSS 39 


II 


K8 




2737 ± 322 


678 ± 10 


300 ± 10 


0.04 


2.07 ±0.32 


1 


ROXs 9D 


III 










< 35 


< 0.007 




3 


WSB 38, VSS 27 


III 


G4 








< 20 


< 0.004 




1 


DoAr 28, HBC 261 


II 


K5 








< 75 


< 0.01 




1 


WL 2, GY 128 


FS 










30 ±10 


0.006 




4 


WL 18, GY 129 


II 


K7 








85 ±10 


0.02 




4 


IRS 23, WL 21 


III 










< 30 


< 0.006 




4 


CRBR 51 


FS 










70 ±15 


0.01 




4 


SR 24, HBC 262 


FS 


Kl 


1271 ± 32 


1263 ± 282 


530 ±9 


230 ± 11 


0.03 


1.29 ±0.23 


3 


WL 14, GY 172 


II 


M4 








30± 10 


0.006 




4 


IRS 27, WL 22 


1 






< 438 


70 ±5 


40 ±10 


0.004 


1.99 ±0.96 


4 


WL 16, YLW 5 


FS 










< 6 


< 0.001 




4 


IRS 29, WL 1 


FS 






< 702 


15 ±4 


< 20 


0.0009 


< 6.05 


4 


LFAM 26, GY 197 


I 










100 ± 10 


0.02 




4 


WL 17, GY 205 


I 






< 867 


100 ± 7 


75 ± 10 


0.006 


0.84 ± 0.65 


4 


WL 10, GY 211 


II 


K8 








< 60 


< 0.01 




4 


GY 213 


II 










40± 10 


0.008 




4 


EL 29, WL 15 


I 






2024 ± 109 


408 ±6 


95 ±10 


0.03 


2.82 ± 0.32 


4 


SR 21, EL 30 


II 


G3 


2761 ± 57 


1896 ± 268 


397 ±6 


95 ±15 


0.03 


2.27 ±0.22 


4 


IRS 31, WL 9 


II 










< 10 


< 0.002 




4 


GY 224 


I 










50± 10 


0.01 




4 


WL 19, GY 227 


II 










30 ± 10 


0.006 




4 


WL 11, GY 229 


II 










< 10 


< 0.002 




4 


WSB 45, HBC 640 


III 


M5 








< 20 


< 0.004 




1 


IRS 33, GY 236 


II 










105 ± 20 


0.02 




4 


WSB 46, HBC 641 


II 


M2 








< 20 


< 0.004 




1 


IRS 34, GY 239 


II 










< 20 


< 0.004 




4 


IRS 36, GY 241 


FS 










< 15 


< 0.003 




4 


WL 20, GY 240 


I 






400 ± 112 


182 ± 6 


47 ±1 


0.01 


1.39 ± 0.40 


5 


IRS 37, GY 244 


II 


M4 




< 591 


93 ±8 


< 10 


0.006 


1.03 ±0.69 


4 


IRS 38, WL 5 


III 


F7 




< 522 


80 ±4 


< 10 


0.005 


< 4.89 


4 


GY 245 


FS 










60 ±5 


0.01 




4 
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Table 1 — Continued 









Fi,(350/im) 


F„(450 /um) 


F^(850 /um) 


F^(1.3mm) 


Md 






Object 


SED 


SpT 


[mjy] 


[mJy] 


[mJy] 


[mJy] 


[Me] 


a 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


TO C on 11 TT A 

IHb oy, WL 4 


II 


M2 




< 480 


63 ± 5 


< 15 


0.01 


1.75 it 0.9o 


4 


TT» O At IUTT O 

1Kb 41, WL 3 


II 






< 582 


115 ± 4 


< 60 


0.007 


< 2.55 


4 


bri 12, HBC zbo 


III 


MO 






< 26 


< 20 


< 0.001 




2,4 


IKo 42, LtY zoz 


II 










< 60 


< 0.01 




4 


WL b, GY 254 


FS 










< 20 


< 0.004 




4 


VBDljr 22, K(JXs 2d 


III 


K6 








< 30 


< 0.006 




1 


WbB 49 


II 


M4 








< 25 


< 0.005 




1 


L-KBK So 


II 










150 it 10 


0.03 




4 


■\/'T 117" 1 1;/^ oco 

YLW loO, GY 2b2 


II 


Ml 








60 it 5 


0.03 




4 


IHb 43, YLW 15 


I 






1811 ± 358 


246 ± 7 


80 it 10 


0.02 


2.99 ± 0.38 


4 


hjL ol, WL la 


II 


MO 








< 10 


< 0.002 




4 


IHo 44, LrY zby 


T 
1 










60 it 10 


0.01 




4 


TTiT QO TOO Ad 

hjL 32, mo 45 


II 


K7 








< 50 


< 0.01 




4 


TOO A£i /^■\/ OT /I 

1Kb 4d, GY 274 


T?0 

I* b 










45 it 10 


0.009 




4 


T7*T OO TOO /I '7 

L/L 3o, 1Kb 47 


II 










< 20 


< 0.004 




4 


CjY 284 


II 


M3 








130 it 10 


0.02 




4 


KUAs 25, CjY 292 


II 


M2 








30 it 5 


0.007 




4 


KUAS oUA 


TTT 
ill 










< 10 


< 0.002 




3 


TOO (4 (^"V or\A 

IHb 48, (jY 3U4 


I 






< 1566 


180 ± 9 


60 it 10 


0.01 


1.91 ± U.Dl 


4 


TOO (^"V QAC 

IKb 5U, CjY 3Ub 


III 


M4 








< 10 


< 0.002 




4 


TOO /in r^v or\Q 
1Kb 49, CjY oUo 


TT 
11 


1\8 




< 1041 


52 it 5 


25 i 5 


0.003 


1.71 di 1.10 


1 


DoAr 32, WbB 51 


II 


K6 








< 45 


< 0.009 




1 


T~i^ A OO 11 rO D CO 

DoAr 33, WbB 53 


II 


K4 




< 678 


79 ± 7 


40 it 10 


0.005 


1.60 ± 0.82 


1 


11 TOO CO /^"V/" O-l A 

WbB 52, CjY 314 


II 


K5 








51 it 10 


0.007 




6 


TOO CI /~*'\7' Q1 C 

IKb 51, CjY 315 


II 






< 678 


171 ± 6 


110 it 10 


0.01 


< 2.17 


4 


OD n AO or\^ 

bK y. Ab 21)7 


II 


K5 






< 25 


15 it 5 


0.002 


< 1.20 


1 


L/L 3d, VbbG 14 


II 


A 'y 
A7 








< 10 


< 0.002 




3 


TO o c /I o^o 

1Kb 54, KjY 378 


FS 










30 it 10 


0.006 




4 




TTT 
ill 


T^8 








^ oc; 

^ 20 


<, U.UUi 




1 o 
1,2 


SR 10, HBC 265 


II 


M2 








< 25 


< 0.005 




1 


WSB 60, YLW 58 


II 


M4 


832 ± 67 


651 ± 161 


149 ±7 


89 it 2 


0.03 


1.59 it 0.24 


5 


SR 20, HBC 643 


HI 


G7 






< 25 


< 35 


< 0.002 




1 


SR 13, HBC 266 


II 


M4 




740 ± 200 


118 ±6 


60 it 10 


0.01 


2.30 it 0.42 


1 


WSB 63 


II 


M2 








< 25 


< 0.005 




1 


ROXs 35A 


HI 


K3 








< 29 


< 0.006 




3 


VSS 44, ROXs 35B 


HI 


G5 








< 19 


< 0.004 




3 


ROXs 39 


HI 


K6 








< 30 


< 0.006 




1 


DoAr 42, HBC 267 


HI 


MO 








< 30 


< 0.006 




1 


ROXs 42B 


HI 


MO 








< 45 


< 0.009 




1 


ROXs 42C 


II 


K6 








< 30 


< 0.006 




1 


ROXs 43A 


II 


GO 








< 35 


< 0.007 




1 


ROXs 43B 


HI 


K5 








< 17 


< 0.003 




3 


IRS 60 


II 










< 50 


< 0.01 




1 


DoAr 44, HBC 268 


II 


K3 




< 414 


181 ±6 


105 it 11 


0.01 


1.28 it 0.59 


3 



Table 1 — Continued 









F^(350/im) 


F^(450 Atm) 


F^(850 Atm) 


-Fi,(1.3 mm) 


Md 






Object 


SED 


SpT 


[mjy] 


[mJy] 


[mJy] 


[mJy] 


[Mq] 


a 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IRS 63, L1709 B 


I 




6703 ± 162 


4729 ± 199 


1017 ± 7 


370 ± 10 


0.06 


2.28 ±0.33 


1 


ROXs 45B 


HI 










< 19 


< 0.004 




3 


ROXs 45C 


HI 


K5 








< 19 


< 0.004 




3 


ROXs 45D 


HI 


KO 








< 19 


< 0.004 




3 


IRS 67, L1689 S 


I 






1687 ± 207 


394 ±9 


150 ± 10 


0.02 


2.28 ±0.33 


1 


DoAr 51, HBC 647 


HI 


MO 








< 20 


< 0.004 




1 


IRS 69 


II 










< 75 


< 0.01 




1 


DoAr 52, HBC 648 


II 


M2 








< 55 


< 0.01 




1 


RNO 90, HBC 649 


II 


G5 




557 ± 83 


162 ± 4 


25 ±5 


0.005 


2.40 ±0.43 


1 


RNO 91, HBC 650 


FS 


Ml 




1524 ± 87 


365 ± 9 


90 ±10 


0.02 


2.64 ±0.32 


1 


EL 49 


II 










58 ± 5 


0.01 




3 


Wa Oph 4, HBC 652 


II 


K4 








< 13 


< 0.003 




1 


Wa Oph 6, HBC 653 


II 


K6 




979 ± 297 


379 ± 7 


130 ± 10 


0.08 


2.07 ±0.40 


1 


Wa Oph 5, HBC 654 


II 


M2 








< 25 


< 0.005 




1 


AS 209, HBC 270 


II 


K5 


3000 ± 72 


1915 ± 206 


551 ± 10 


300 ± 10 


0.05 


1.80 ±0.22 


1 



Note. — Col. (1): Object name(s). Col. (2): SED classification type (FS = Flat Spectrum); see §3.3 for details. Col. (3): Stellar 
spectral type from the literature (see §4.2). Col. (4): 350 fim flux density. Col. (5): 450 fim flux density. Col. (6): 850 fim flux density. 
Col. (7): 1.3 mm flux density. Col. (8): Disk mass (see §3.1). Col. (9): Submillimeter continuum s lope (see §3.2). Col. (10): Notes 
on individual sources as follows: (1) The 1.3 mm flux densities are from I An dre &: Montmerl3 lll994h . The la rms uncertainties were 
typically assumed to be ~10mjy. (2) 80 fim flux densities from Ijensen et al. I l ll996h are listed in Col. (6). (3) The 1.3mm flux 
densities are from lNiirnberger et all i ll998h . (4) The 1.3 mm flux densities are from lMotte et al. I lll998l ). (5) The 1.3 mm flux densities 
for WL 20 (S) and WSB 60 are from the interferometric survey of [Andrews fc Williamsl ||2007| ). (6) The 1.3 mm flux densities are 
from lstanke et"aL 1 1|2006| ). All flux densities are measured in units of mJy. Upper limits are taken at the 3 a confidence level. Quoted 
errors are the 1 cr rms noise levels and do not include systematic errors in the absolute flux calibration (~25% at 350 and 450 fim, 
~10% at 850 /xm, and ~20% at 1.3 mm). 
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Table 2. Results of SED Fits 





Ti 














Ti 




Ma 








Object 


[K] 


g 


[10-3 Mq] 




dof 


Notes 


Object 


[K] 


1 


[10-3 Mq] 




dof 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


AS 205 


304 


0.55 


33 ± 5 


2.6 


7 


1,2,3 


IRS 49 


142 


0.58 


2.9 ±0.6 


0.5 


3 


1,2,5 


SR 4 


181 


0.47 


3.8 ±0.5 


0.6 


4 


1,2,4 


WSB 52 


162 


0.52 


7±3 


0.1 


3 


1,2,5,7 


GSS 26 


157 


0.63 


45 ± 11 


3.8 


4 


1,2,4 


SR 9 


150 


0.56 


1.9 ±0.9 


8.7 


2 


1,2,5 


EL 20 


157 


0.69 


33 ± 9 


3.0 


3 


1,2,5 


WSB 60 


99 


0.57 


32 ±8 


4.0 


6 


1,2,6 


DoAr 24E 


218 


0.63 


8± 1 


0.3 


5 


1,2,5,6 


SR 13 


189 


0.63 


14 ±3 


1.0 


6 


1,2,3 


EL 24 


229 


0.65 


162 ± 62 


1.5 


8 


1,3,5,6 


DoAr 44 


169 


0.54 


10 ±2 


1.2 


3 


1,2,3 


EL 26 


116 


0.70 


5 ± 3 


2.2 


2 


1,2,5 


RNO 90 


215 


0.59 


4.7 ±0.7 


5.7 


5 


1,2,3 


IRS 39 


104 


0.69 


14 ±5 


2.0 


4 


1,2,3,6 


Wa Oph 6 


157 


0.65 


81 ±23 


1.5 


5 


3,8 


YLW 16C 


120 


0.63 


26 ± 10 


0.2 


2 


1,2,5 


AS 209 


230 


0.63 


52± 15 


1.0 


4 


3 


GY 292 


143 


0.63 


7±3 


2.3 


2 


1,2,6 
















Note. — 


Col. (1): 


Object 


name. Col. (2): 


Best-fit value of the temperat 


ure at 


1 AU 


in K (typical 


errors 


are ± 


a few K). 



Col. (3): Best-fit value of the radial power-law index of the temperature profile (typical errors are ±0.02). Col. (4): Best-fit 
value and uncertainty of the disk mass. Col. (5): Reduced statistic. Col. (6): Number of degrees of freedom in the fit 
(i.e., number of data points minus number of fitted parameters [=3]). Col. (7): Notes on literature sources for the SED data as 
follows: (1) Spitzer Space Telescope IRAC flux densities (3.6, 4.5, 5.8, 8 fim) from L. Allen (2007, private communication). (2) 
Spitzer Space T elescope MI PS fiux density (24 fivci) from L. Allen (2007, private comm unication) . (3) IRAS flux densities (12, 
25, 60 /.tm) from lWeaver fc Jo ne s (199g). (4) Mid- in frared photometry (10, 20 lira) from lwilking et al. (1989). (5) Mid-in frared 
photometry (10, 12 /im) from lBarsoii7et al. | |2005l'). (6) ISO p hotometry (6.7 and 14.3 um) from iBontemps et al.l ll200ll'). (7) 
Mid-in frared photometry (lO. 20 um") from lcreene et all l ll994l ). (8) Spitzer IRAC and MIPS photometrv from [Padgett et al.l 
ll2006h . 
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Table 3. Statistical Comparison of p Oph and Taurus Submillimeter Properties 



Samples 


[%] 


P{a) 
[%] 


Regional Comparison 






Tau vs. Oph (all) 


35-47 


9-62 


Tau I vs. Oph I 


20 - 72 


54 - 55 


Tau II vs. Oph II 


7-10 


3-65 


Tau 111 vs. 0])li ill 


9-34 




EvoluUonary Comparison 
Oph I vs. Oph II 
Oph II vs. Oph III 
Tau I vs. Tau II 
Tau II vs. Tau III 


> 99.75 

> 99.9999 

> 99.97 

> 99.9999 


92-95 
> 99.91 



Note. — Results of survival analysis two-sample statisti- 
cal tests comparing various subsamples of the submillimeter 
properties of Taurus (Paper I) and p Oph YSOs. The listed 
values P correspond to the probabilities that the two sub- 
samples are drawn from different parent distributions. The 
probability ranges are from different statistical tests. 



-27- 



Table 4. Median Submillimeter Properties of YSOs 





(Ma) 


{Ma /Ah) 




Sample 


[Mq] 


[%] 


(a) 


Class I 


0.015 


3.8 


2.3 


Class II 


0.005 


0.9 


1.9 


transition^ 


0.0009 


0.3 




Class III 


< 5 X 10-^ 







^Transition objects are defined here to be 
Class III sources (i.e., n > 1.6) that have 
submillimeter detections. There is 1 such ob- 
ject in p Oph (WL 5) and 7 others in Taurus 
(CIDA-8, CoKu Tau/4, FW Tau, FY Tau, HQ 
Tau, LkHo 332/Gl, and V807 Tau). The sub- 
milli meter detecti o n of WL 5 may be anoma- 
lous: iMotte et al. I il998l) showed that the orie- 
inal quoted 1.3 mm detecti on from earlier work 
iAndre fc Montmerlj Il994) was actually back- 
ground emission. The Class III sample median 
is based on the stacked average limit derived 
in 54.1. 
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Fig. 1. — Histograms of some relevant properties for the p Oph sample obie cts. The left and middle 



panels show stellar masses and ages, inferred as described in §4.2 using the lD'Antona &: Mazzitelh 



(jl997l ) theoretical pre-main— sequence evolution models. The right panel shows the infrared SED 
index (defined as n where vFjj oc z^"; see §3.3) that is commonly used as an evolutionary diagnostic 
for circumstellar dust. YSOs with an index less than —0.3 are Class I sources, between —0.3 and 
1.6 are Class II sour ces, and > 1.6 are C lass III sources (SED classification borders are shown as 
dotted vertical lines: iGreene et al.lll994l ). 
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Md [Mq] Ma [Mo] 



Fig. 2. — Correlations between the 850 //m {left) or 1.3 mm [right) flux densities and the best-fit 
disk masses inferred from fitting SEDs with a simple disk structure model. The results for p Oph 
{circles: see Table [2]) and Taurus {asterisks: see Paper I) YSOs are shown together (the latter flux 
densities have been scaled to c? = 150 pc for ease of comparison). Overlaid as a solid line in each 
panel is the generic relationship described by Eqn. 2 for = 20 K, used to calculate disk masses 
(or upper limits) for the remainder of the sample. 
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Fig. 3. — The cumulative distribution of disk masses for the p Oph sources Usted in Table [H 
constructed using the Kaplan-Meier estimator to include non-detections. The ordinate represents 
the probability that a sample disk has a mass greater than the abscissa value. The median disk 
mass for the full sample is ~0.002 Mq. 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 
a 

Fig. 4. — The cumulative distribution of submillimeter colors {a, defined as Fy cc u^) for the p Oph 
sources listed in Table [TJ constructed using the Kaplan- Meier estimator to include non-detections. 
The ordinate represents the probability that a sample disk has a spectral index greater than the 
abscissa value. The median value for the full sample is (a) 2.0. 
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Fig. 5. — Direct comparisons of submillimeter properties in the p Oph (black lines) and Taurus 
(grayscale) regions. The panels on the left and right show the cumulative distributions of and 
a, respectively, for various subgroups of YSOs in each region. The top panels show the full samples 
described here and in Paper I. The subsequent panels below show the distributions for Class I, 
II, and III sources. The distributions of these submillimeter properties in the p Oph and Taurus 
regions are essentially indistinguishable. However, there are statistically significant changes in the 
Md and a distributions along the Class I ^ II — > III evolution sequence. 
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Fig. 6. — Individual disk masses (top) and submillinictcr colors (a: bottom) as a function of the 
infrared SED index measured between ~2— 25 ;um (ra) for the combined p Oph (circles) and Taurus 
(asterisks) samples. Open triangles represent upper limits. Both panels show correlations such that 
Md and a decrease along the empirical evolution sequence (as n increases). 
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MJM. 

Fig. 7. — The cumulative distributions of the disk-to-star mass ratio for Class II sources in p Oph 
{lines) and Taurus {grayscale). While the p Oph sources appear to generally have slightly larger 
mass ratios than their Taurus counterparts, two-sample tests on the censored data indicate that 
the difference is not statistically significant. The median ratios are ~1.0% and 0.8% for p Oph and 
Taurus, respectively. 
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Fig. 8. — Disk masses, disk-to-star mass ratios, and submillimeter colors as a function of the stellar 
age for the combined p Oph {circles) and Taurus {asterisks) samples. Open triangles are upper 
limits. Correlation tests on the censored data show that the disk-to-star mass ratio decreases 
significantly with time, although similar trends for and a are not clear. 
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Fig. 9. — Schematic diagram showing the evolution of the median disk masses across the empirical 
sequence defined by the shape of the infrared SED (see also Table H]). The transition objects 
("trans.") are defined as Class III sources with firm submillimeter detections, and the Class III 
limit corresponds to the stacking analysis value determined in §4.1. Dashed horizontal lines mark 
standard reference masses. 
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Fig. 10. — The differential distributions of and a for Class II sources that were firmly detected 
at submillimeter wavelengths in the combined p Oph and Taurus surveys. Uncertainties correspond 
to the counting errors in each bin. These distributions are representative of the known properties 
of disks that may represent the initial conditions for the planet formation process. 
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Fig. 11. — A comparison of disk masses inferred in two independent ways for p Oph (circles) and 
Taurus (asterisks) Class II disks. The abscissa uses gas accretion rates from optical/infrared spectra 



and some assumptions about the accretion history (Md ~ 2t^M: lHartmann et al.lll998l : lNatta et al 



20061 ) and the ordinate the dust masses described in §3.1 and Paper I. The solid line marks the 



expected behavior if the values were equal. The dotted line bisects the sample, and indicates that 
dust masses underestimate the accretion masses by nearly an order of magnitude. Particle growth 
up to meter size scales may explain the discrepancy, due to a significant over-estimate of the disk 
opacity. 



